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Abstract 

The purpose of this work is to understand the zero temperature phases, and 
the phase transitions, of Heisenberg spin systems which can have an extensive, 
spontaneous magnetic moment; this entails a study of quantum transitions 
with an order parameter which is also a non-abelian conserved charge. To 
this end, we introduce and study a new class of lattice models of quantum 
rotors. We compute their mean-field phase diagrams, and present continuum, 
quantum field-theoretic descriptions of their low energy properties in different 
regimes. We argue that, in spatial dimension d = 1, the phase transitions in 
itinerant Fermi systems are in the same universality class as the corresponding 
transitions in certain rotor models. We discuss implications of our results for 
itinerant fermions systems in higher d, and for other physical systems. 
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I. INTRODUCTION 



Despite the great deal of attention lavished recently on magnetic quantum critical phe- 
nomena, relatively little work has been done on systems in which one of the phases has 
an extensive, spatially averaged, magnetic moment. In fact, the simple Stoner mean- field 
theory of the zero temperature transition from an unpolarized Fermi liquid to a ferro- 
magnetic phase is an example of such a study, and is probably also the earliest theory of 
a quantum phase transition in any system. What makes such phases, and the transitions 
between them, interesting is that the order parameter is also a conserved charge; in sys- 
tems with a Heisenberg 0(3) symmetry this is expected to lead to strong constraints on the 
critical field theories [0]. As we will discuss briefly below, a number of recent experiments 
have studied systems in which the quantum fluctuations of a ferromagnetic order parameter 
appear to play a central role. This emphasizes the need for a more complete theoretical 
understanding of quantum transitions into such phases. In this paper, we will introduce 
what we believe are the simplest theoretical models which display phases and phase transi- 
tions with these properties. The degrees of freedom of these models are purely bosonic and 
consist of quantum rotors on the sites of a lattice. We will also present a fairly complete 
theory of the universal properties of the phases and phase transitions in these models, at 
least in spatial dimensions d > 1. Our quantum rotor models completely neglect charged 
and fermionic excitations and can therefore probably be applied directly only to insulating 
ferromagnets. However, in = 1, we will argue that the critical behaviors of transitions in 
metallic, fermionic systems are identical to those of the corresponding transitions in certain 
quantum rotor models. 

We now describe the theoretical and experimental motivation behind our work: 
(i) The Stoner mean field theory of ferromagnetism in electronic systems in fact con- 
tains two transitions: one from an unpolarized Fermi liquid to a partially polarized itinerant 
ferromagnet (which has received some recent experimental attention 0), and the second 
from the partially polarized to the saturated ferromagnet. A theory of fluctuations near the 
first critical point has been proposed but many basic questions remain unanswered , 
especially on the ordered side |^ (there is no proposed theory for the second transition, 
although we will outline one in this paper). It seems useful to examine some these issues 
in the simpler context of insulating ferromagnets. Indeed, as we have noted, we shall argue 
below that in ci = 1, certain insulating and itinerant systems have phase transitions that are 
in the same universality class. 



[ii) Many of the phases we expect to find in our model also exist in experimental com- 
pounds that realize the so-called "singlet-triplet" model 0] . These compounds were studied 
many years ago with a primary focus on finite temperature, classical phase transitions; we 
hope that our study will stimulate a re-examination of these systems to search for quantum 
phase transitions 

(iii) All of the phases expected in our model (phases A-D in Section |iB| below), occur in 
the Lai_xSrxMnO^ compounds These, and related compounds, have seen a great deal 
of recent interest for their technologically important "colossal magnetoresistance" . 



[iv) Recent NMR experiments by Barrett et. al. |Tl[] have studied the magnetization of a 



quantum Hall system as a function of both filling factor, i/, and temperature, T, near v = 1. 
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The T = state at z/ = 1 is a fully polarized ferromagnet [p!^ -[T5[] and its finite temperature 
properties have been studied from a field-theoretic point of view |T^. More interesting for 
our purposes here is the physics away from z/ = 1: Brey et. al. have proposed a variational 
ground state consisting of a crystal of "skyrmions"; this state has magnetic order that is 
canted |I8[, i.e. in addition to a ferromagnet moment, the system has magnetic order (with 
a vanishing spatially averaged moment) in the plane perpendicular to the average moment. 
A phase with just this structure will appear in our analysis, along with a quantum-critical 
point between a ferromagnetic and a canted phase. Although our microscopic models are 
quite different from those appropriate for the quantum Hall system, we expect the insights 
and possibly some universal features of our results to be applicable to the latter. 

There has also been some interesting recent work on the effects of randomness on itinerant 
ferromagnets . This paper shall focus exclusively on clean ferromagnets, and the study of 
the effects of randomness on the models of this paper remains an interesting open problem; 
we shall make a few remarks on this in Section |V11 3| 

In the following subsection we will introduce one of the models studied in this paper, 
followed by a brief description of its phases in Section [TBI . Section | will conclude with an 
outline of the remainder of the paper. 



A. The Model 

We introduce the quantum rotor model which shall be the main focus of the paper; 
extensions to related models will be considered later in the body of the paper. On each site 
z of a regular lattice in d dimensions there is a rotor whose configuration space is the surface 
of a sphere, described by the 3-component unit vector hi^ {fi = 1,2,3 and J^fifi'l^ = 1); the 
caret denotes that it is a quantum operator. The canonically conjugate angular momenta 
are the Lj^, and these degrees of freedom obey the commutation relations (dropping the site 
index as all operators at different sites commute) 

[n^,n^]=0 ; [Lf,, L^] = ie^^xLx ; [Lf„hu] = ie^^xfix (1.1) 

As an operator on wavefunctions in the configuration space, is given by 

d 

Lf, = -ie^^yxnuT, — (1-2) 
dnx 

We will be interested primarily in the properties of the Hamiltonian 

^ = f II + « J ~ {J^it^^m + KLi^Lj^ + M {fii^Ljf, + nj^U^)) (1.3) 



<*j> 



where there is an implied summation over repeated /x indices, and < ij > is the sum over 
nearest neighbors, and the couplings g, a, J, K, M are all positive. All previous analyses 
of quantum rotor models |]20|- p^] have focussed exclusively on the case K = M = 0, and the 



novelty of our results arises primarily from nonzero values of the new K, M couplings. A 
crucial property of H is that the 3 charges 
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Q,=T.L^, (1.4) 
i 

commute with it, and are therefore conserved. Indeed, H is the most general Hamiltonian 
with bihnear, nearest neighbor couphngs between the and operators, consistent with 
conservation of the Q^. We have also included a single quartic term, with coefficient ga, 
but its role is merely to suppress the contributions of unimportant high energy states. 

Discrete symmetries of Ti. will also be important in our considerations. Time-reversal 
symmetry, T is realized by the transformations 

T : L^-^ -Lf, -fif, (1.5) 

Notice that the commutators ( |1 . 1| ) change sign under T, consistent with it being an anti- 
unitary transformation. All the models considered in this paper will have T as a symmetry. 
For the special case M = 0, we also have the additional inversion symmetry V: 

V : L^-^ -hf, (1.6) 

The presence of V will make the properties of the M = system somewhat different from 
the M 7^ case. We will see later that V is related to a discrete spatial symmetry of the 
underlying spin system that H models. 

The utility of H does not lie in the possibility of finding an experimental system which 
may be explicitly modeled by it. Rather, we will find that it provides a particularly simple 
and appeahng description of quantum phases and phase transitions with a conserved order 
parameter in a system with a non-abelian symmetry. Further, we will focus primarily on 
universal properties of H, which are dependent only on global symmetries of the states; 
these properties are expected to be quite general and should apply also to other models 
with the same symmetries, including those containing ordinary Heisenberg spins. 

To help the reader develop some intuition on the possibly unfamiliar degrees of freedom in 
H, we consider in Appendix ^ a general double-layer Heisenberg spin model containing 



both inter- and intra- layer exchange interactions. We show that, under suitable conditions, 
there is a fairly explicit mapping of the double layer model to the quantum rotor Hamiltonian 
H. Under this mapping we find that each pair of adjacent spins on the two layers behaves 
like a single quantum rotor. In particular ~ S'^^ + Si,^ and ~ S'a^ — 5*6^ where a, b are 
the two layers and the are Heisenberg spins. Notice also that P is a layer-interchange 
symmetry. 

Before turning to a description of the ground state of H, it is useful to draw a parallel to 
another model which has seen a great deal of recent interest — the boson Hubbard model . 



The latter model has a single conserved charge. A*";, the total boson number, associated with 
an abelian global U{1) symmetry. In contrast, the quantum rotor model H has the 3 charges 
Qfj^, and a non-abelian global 0(3) symmetry. As we will see, the non-abelian symmetry 
plays a key role and is primarily responsible for the significant differences between H and 
the boson Hubbard model. It is also useful to discuss a term-by-term mapping between 
H and the boson Hubbard model. The terms proportional to g in H are analogous to the 
on-site Hubbard repulsion in the boson model. The latter model also has an on-site chemical 
potential term which couples linearly to Nb, but such a term is prohibited by symmetry in 
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the non-abelian rotor model. The J term in H has an effect similar to the boson hopping 
term, while the K term is like a nearest-neighbor boson density- density interaction. There 
is no analog of the M term in the boson Hubbard model. 

B. Zero temperature phases of H 

We show in Fig |I| the zero temperature (T) phase diagram of H in the J plane at fixed 
(7, a, and M. This phase diagram was obtained using a mean-field theory which becomes 
in exact in the limit of large spatial dimensionality {<£) ; however, the topology and general 
features are expected to be valid for all c? > 1. The d = 1 case will be discussed separately 
later in the paper; in the following discussion we will assume d > 1. We will also assume 
below that V symmetry is absent, unless otherwise noted. Throughout this paper we will 
restrict consideration to parameters for which the ground state of H are translationally 
invariant ground states — this will require that M not be too large. 

There are four distinct classes of phases: 

(A) Quantum Paramagnet: 

This is a featureless spin singlet and there is a gap to all excitations. The 0(3) symmetry 
remains unbroken, as 

{h,) = (l,) = 0. (1.7) 

Clearly this phase will always occur when g is much bigger than all the other couplings. 

(B) Quantized Ferromagnets: 

These are ordinary ferromagnets in which the total moment of the ground state is quan- 
tized in integer multiples of the number of quantum rotors (extensions of H in which the 
quantization is in half-integral multiples will be considered later in this paper). The fer- 
romagnetic order parameter chooses a direction in spin space (say, z), but the symmetry 
of rotations about this direction remains unbroken. The ground state therefore has the 
expectation values 

(l,) = integer ^ ; (n,) ^ (1.8) 

The value of {"hz) is not quantized and varies continuously as K, J, and M are varied; for the 
system with V symmetry (M = 0) we will have (fiz) = in this phase. If we consider each 
quantum rotor as an effective degree of freedom representing a set of underlying Heisenberg 
spins (as in the double-layer model of Appendix 0), then (n^) determines the manner in 
which the quantized moment is distributed among the constituent spins. The low-lying 
excitation of these phases is a gapless, spin-wave mode whose frequency a; ~ fc^ where k is 
the wavevector of the excitation. 

That these phases occur is seen as follows: First consider the line J = 0, M = 0. For very 
small K, it is clear that the ground state is a quantum paramagnet. As K is increased, it is 
easy to convince oneself (by an explicit calculation) that a series of quantized ferromagnet 
phases with increasing values of (^Lz'^ = integer get stabilized. Further in this simple limit 
the exact ground state in each one of these phases is just a state in which each site is put 
in the same eigenstate of and L^. There is a finite energy cost to change the value of 
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at any site. Now consider moving away from this hmit by introducing small non-zero 
values of J and M. These terms vanish in the subspace of states with a constant value of 
Lf so we need to consider excitations to states which involve changing the value of Lf at 
some site. As mentioned above such states are separated from the ground state by a gap. 
Consequently, though the new ground state is no longer the same as at J = 0, M = 0, it's 
quantum numbers , in particular the value of Qz = J2iLiz, are unchanged. The stability of 
the quantized ferromagnet phases up to finite values of J and M implies the existence of 
direct transitions between them which is naturally first order. In general a non-zero value 
of M will also lead to a non-zero value of (n^). 

All of this should remind the reader of the Mott-insulating phases of the boson Hubbard 
model [^. In the latter, the boson number, n^, is quantized in integers, which is the analog 



of the conserved angular momentum of the present model. However, the properties of 
the Mott phases are quite similar for all values of {fib), including {fih) = 0. In contrast, for 
the 0(3) rotor model, the case (^Lz'^ = (the quantum paramagnet, which has no broken 

symmetry and a gap to all excitations) is quite different from (^Lz'^ = integer 7^ cases 
(the quantized ferromagnets, which have a broken symmetry and associated gapless spin- 
wave excitations). As we argued above and as shown in Fig the lobes of the quantized 
ferromagnets and the quantum paramagnet are separated by first-order transitions, while 



there were no such transitions in the phase diagram for the boson Hubbard model in Ref p5 



this is an artifact of the absence of off-site boson attraction terms analogous to the K term 



in Ref ||25|, and is not an intrinsic difference between the abelian and non-abelian cases. 
(C) Neel Ordered Phase: 
In this phase we have 



(n^) ^ 0, while {L^) = 0. (1.9) 

This occurs when J is much bigger than all other couplings. We refer to it as the Neel 
phase because the spins in the bilayer model are oriented in opposite directions in the two 
layers. More generally, this represents any phase in which the spin-ordering is defined by a 
single vector field (n^) and which has no net ferromagnetic moment. There are 2 low-lying 
spin- wave modes, but they now have a linear dispersion u; ~ [|^. The order parameter 
condensate (n^) does not have a quantized value and varies continuously with changes in 
the couplings. 

From the perspective of classical statistical mechanics, the existence of this Neel phase is 
rather surprising. Note that there is a linear coupling MLj^rij^ in H between the ri^ and 
fields. In a classical system, such a coupling would imply that a non-zero condensate of 
must necessarily accompany any condensate in n^. The presence of a phase here, in which 
= despite (ri^) 7^ 0, is a consequence of the quantum mechanics of the conserved 
field L^. 

(D) Canted Phase: 

Now both fields have a condensate: 

(n^) ^ and (4) ^ 0. (1.10) 

The magnitudes of, and relative angle between, the condensates can all take arbitrarily 
values, which vary continuously as the couplings change; however if P is a symmetry (if 
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M = 0) then (n^) is always orthogonal to \L^l/■ The 0(3) symmetry of H is completely 
broken as it takes two vectors to specify the orientation of the condensate. In the bilayer 
model, the spins in the two layers are oriented in two non-collinear directions, such that 
there is a net ferromagnetic moment. The low lying excitations of this phase consist of 2 
spin wave modes, one with u; ~ fc, while the other has u; ~ fc^. This phase is the analog 
of the superfluid phase in the boson Hubbard model. In the latter model the conserved 
number density has an arbitrary, continuously varying value and there is long range order 
in the conjugate phase variable; similarly here the conserved (Lf^ has an arbitrary value 
and the conjugate field has a definite orientation. 

At this point, it is useful to observe a parallel between the phases of the quantum rotor 
model and those of a ferromagnet Fermi liquid. This parallel is rather crude for c/ > 1, but, 
as we will see in Sees and |VI1] , it can be made fairly explicit in c? = 1 where we believe 
that the universality classes of the transitions in the rotor model and the Fermi liquid are 
identical. The quantized ferromagnetic phases B are the analogs of the fully polarized phase 
of the Fermi liquid in which the spins of all electrons are parallel. The canted phase D has 
a continuously varying ground state polarization, as in a partially polarized Fermi liquid. 
Finally the Neel phase C is similar to an unpolarized Fermi liquid in that both have no net 
magnetization and exhibit gapless spin excitations. 

Finally, we point out an interesting relationship between phases with a non-quantized 
value of the average magnetization (as in phase D of the rotor model), and phases with 
vanishing average magnetization (phases A and C). Notice that, in both the Fermi liquid 
and the quantum rotor model, it is possible to have a continuous transition between such 
phases. However, in both cases, such a transition only occurs when the phase with vanishing 
magnetization has gapless spin excitations; in phase C of the rotor model we have the gapless 
spin waves associated with the Neel long range order, while in the unpolarized Fermi liquid 
we have gapless spin carrying fermionic quasiparticles. Phase A of the rotor model has 
no gapless excitations, but notice that there is no continuous transition between it and 
the partially polarized canted phase D. We believe this feature to be a general principle: 
continuous zero temperature transitions in which there is an onset in the mean value of a 
non-ahelian conserved charge only occur from phases which have gapless excitations. 

The outline of the remainder of the paper is as follows. We will begin in Section |I| by 
discussing the mean field theory which produced the phases described above. In Section ^ 
we will perform a small fluctuation analysis of the low-lying excitations of the phases. We 
will turn our attention to the quantum phase transitions in (i > 1 in Section ^V|, focusing 
mainly on the transition between phases C and D in Section |IV A| , and that between phases 
B and D in Section |1V B| . In Section |V| we shall consider an extension of the basic rotor model 
( |1.3| ): each rotor will have a 'magnetic monopole' at the origin of n space, which causes the 
angular momentum of each rotor to always be non-zero. We will turn our attention to the 
important and distinct physics in (i = 1 in Section |VI[ We will conclude in Section |VII| by 
placing our results in the context of earlier work and discuss future directions for research. 
Some ancillary results are in 5 appendices: we note especially Appendix 0, which contains 
new universal scaling functions of the dilute Bose gas, a model which turns out to play a 
central role in our analysis. 
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II. MEAN FIELD THEORY 



In this section we will describe the mean-field theory in which the phase diagram of Fig [I| 
was obtained. The mean-field results become exact in the limit of large spatial dimensional- 
ity, d. Rather than explicitly discussing the structure of the large d limit, we choose instead 
a more physical discussion. We postulate on every site a single-site mean field Hamiltonian 

Hn.f = \{Ll + a {Liy^ - N,n, - h.L, (2.1) 

which is a function of the variational c- number local fields A^^ and h^. These fields are deter- 
mined at T = by minimizing the expectation value of H in the ground state wavefunction 
of Hmf- The mean-field ground state energy of H is 

Emf = Eo-^ {nX - ^ - MZ {n,)^ (4)^ + N, {n,)^ + (L^)^ (2.2) 

where Eq is the ground state energy of the Hmf, all expectation values are in the ground 
state wavefunction of Hmf, and Z is the co-ordination number of the lattice. We now have 
to minimize the value of Emf over variations in /i^ and A^^. This was carried out numerically 
for a characteristic set of values of the coupling constants. Stability required that the cou- 
pling M not be too large. Further details may be found in Appendix Here we describe 
the behavior of the local effective fields h^, A^^ in the various phases: 

(A) Quantum Paramagnet: 

This phase has no net effective fields A"^ = 0, /i^ = 0. 

(B) Quantized Ferromagnets: 

We now have and ^ with all other components zero. The values of A^^; and 

both vary continuously as the parameters are changed. Nevertheless, the value of (^L^'^ re- 
mains pinned at a fixed non-zero integer. This is clearly possible only because commutes 
with H and Hmf, and Lz is therefore a good quantum number. 

(C) Neel Ordered Phase: 

Like phase B, this phase has ^ and hz with all other components zero, and 
the values of A^^ and hz both vary continuously as the parameters are changed. However 
-^^) = 0; it is now quantized at an integer value which happens to be zero. As a result, 
there is no net ferromagnetic moment in this phase. This unusual relationship between an 
order parameter ^-f'/i^, and its conjugate field /i^, is clearly a special property of the interplay 
between quantum mechanics and conservation laws, and cannot exist in classical statistical 
mechanics systems. 

(D) Canted Phase: 

Now both hn and A"^ are non-zero, and take smoothly varying values with no special con- 
straints, as do their conjugate fields and L^. In systems with V a good symmetry 
(M = 0), we have E/. hf^N,, = 0. 
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III. STRUCTURE OF THE PHASES 



The quantum paramagnet A is a featureless singlet phase with all correlations decaying 
exponentially in both space and imaginary time, and a gap to all excitations. The ferro- 
magnet phases B and the Neel phase C are conventional magnetically ordered phases and 
hardly need further comment here. We describe below the long-wavelength, low energy 
quantum hydrodynamics of the canted phase D. We are implicitly assuming here, and in the 
remainder of this section that d > 1. 

We will study the phase D by accessing it from the Neel phase C. We will analyze 
the properties of D for small values of the uniform ferromagnetic moment; this leads to a 
considerable simplification in the analysis, but the form of the results are quite general and 
hold over the entire phase D — in a later section (Section |I V B ) , we will also access phase D 



from one of the quantized ferromagnetic phases B and obtain similar results. 

We will use an imaginary time, Lagrangian based functional-integral point of view. The 
analysis begins by decoupling the inter-site interactions in 7i by the spacetime dependent 
Hubbard Stratonovich fields A'"^(x,r) and h^{x,T); these fields act as dynamic local fields 
similar to those in Hrnf (Eqn ( p.l| )). We can then set up the usual Trotter product decompo- 
sition of the quantum mechanics independently on each site: this yields the following local 
functional integral on each site (we are not displaying the inter-site terms involving the iV^ 
and fields, as these will be considered later): 



Zl = J T^ni,5{nl - 1) exp 




1 fdn Y 
= ^ ( - ie^uxhunxj - rifj^Nf, (3.1) 

We have ignored, for simplicity, the contribution of the quartic a term in Ti. It is not possible 
to evaluate exactly; for time-independent source fields h^, Nfj_, the evaluation of is of 
course equivalent to the numerical diagonalization that was carried out in Section ||. For 
A^^ = 0, however, we can obtain the following simple formula for the ground state energy 
El = -lim^^oo(l//3) ^ogZ^. 

El = Min [gi{i + l)/2 ~ £h] ; £ > , integer. (3.2) 

The minimum is taken over the allowed values of i, and we have again ignored the a term. 
Note that this is a highly non-analytic function of h: these non-analyticities are directly 
responsible for the lobes of the quantized ferromagnetic phases B. In this section we will 
begin by working in the region of parameters in which El is minimized by ^ = (phases A 
and C); in the vicinity of this region it is permissible to expand in powers of h. Inter-site 
effects will then eventually lead to a phase in which there is a net uniform moment (phase 
D) — this moment will not be quantized and there will be appreciable fluctuations in the 
magnetic moment of each site (this is similar to fluctuations in particle number in a boson 
superfiuid phase). For simplicity we will first present the analysis for the case M = 0. Later 
we will indicate the modifications necessary when M ^ 0. 
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A. Model with V symmetry 



Recall that V symmetry is present when M = 
It is convenient to write an effective action functional in terms of the and fields: 
the form of this functional can be guessed by symmetry and the usual Landau arguments: 



Pn^m^exp j d'^x J^^ dT{Ci + €2)^ 




+ -^hl + ^{hl) 



= + +f{nir + f K)(/.^) + ^{n,h,r (3.3) 

We have temporarily modified from a fixed-length to a "soft- spin" field; this is merely for 
convenience in the following discussion and not essential. Here if^ is an external magnetic 
field whose coupling to the fields is determined by gauge-invariance [^. The reason for 
splitting the Lagrangian into pieces £1, £2 will become clear below. 

We can now look for static, spatially uniform, saddle points of C. This gives three 
different types of solutions, corresponding to the phases A, C, and D (the absence of a 
length constraint on is necessary to obtain all three saddle-points at tree level). The 
values of the and fields at these saddle points are identical in form to those discussed 
for these phases in Section y (with the reminder that we have temporarily specialized to 
M = so that P is a good symmetry). Note that there is no saddle point corresponding to 
the B phases: this is clearly a consequence of ignoring the non-analytic behavior of h'm Z^. 

The low-lying excitation spectrum in the A, C, and D phases can now be determined 
by an analysis of Gaussian fluctuations about the saddle points. While simple in principle, 
such an analysis is quite tedious and involved, especially in phase D. We will therefore not 
present it here; we present instead a more elegant approach in which the answer can be 
obtained with minimal effort. 

Recall that an efficient method of obtaining the properties of the Neel phase is to use 
a non-linear sigma model in which the constraint = 1 is imposed. This eliminates high 
energy states from the Hilbert space associated with amplitude fluctuations, but does not 
modify the low energy spectrum. We introduce here an extension to a hybrid sigma model 
which allows also for the existence of the canted phase D. Notice from the mean- field 
solutions in Sec fT| that the onset of the canted phase D is signaled by the appearance of 
an expectation value of /i^ in a direction orthogonal to the mean direction of n^, while 
the component of /i^ parallel to is zero in both the C (Neel) and D phases (note: the 
last restriction on components of /i^ and parallel to each other requires V symmetry 
and M = 0). This suggests that the important fluctuations are the components of 
perpendicular to ra^, while fluctuations which change the dot product n^/i^ are high energy 
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modes. So we define our hybrid sigma model by imposing the two rotationally invariant 
constraints 



n^n^ = 1 ; n^h^ = (3.4) 

With these constraints, the degrees of freedom have been reduced from the original 6 real 
fields h^, to 4. Notice that while amplitude fluctuations of have been eliminated, those 
in the components of orthogonal to have not — this is the reason for the nomenclature 
'hybrid' above. As an immediate consequence of ( |3.4| ), all the terms in £2 either become 
constants or modify couplings in £1, and £1 is the Lagrangian of the hybrid sigma model. 

The Lagrangian Ci and the constraints ( ^.4| ) display, in principle, all three phases A, 
C, and D. We now move well away from phase A, assuming that has a well developed 
expectation value along the z direction. We then parametrize deviations from this state by 
the following parametrization of the fields, which explicitly obeys the constraints 



n 



We have two complex fields ip, 0, corresponding to the 4 degrees of freedom in the hybrid 
sigma model. The field ip represents oscillations of the about its mean value: by definition, 
we will have (ip) = in all phases. The field measures the amplitude of /i^ orthogonal to 
the instantaneous value of n^: we have (0) = in the Neel phase C, while (0) 7^ in the 
canted phase D. We now insert ( |3.5| ) into ( p.3| ) and obtain at if = the partition function 
Z„ for our hybrid sigma model 



£.1 = 



dip ^ 



+ Ki (0*^-0^1 +ir2|V^|2 + ir3|V0|2 + r4| 




+ ^ I V(V^*0 + + ^(^*0 + V^0*)' (3.6) 

where = ri + vi — Ki. The Lagrangian £0-1 (£0-2) contains terms that are quadratic 
(quartic) in the fields. Note that, apart from the Ui term, all coupling constants in £0-2 
are related to those in £0-1: these constraints on the couplings encapsulate the rotational 
invariance of the underlying physics. The field ip, which represents fluctuations of iV^ about 
its average value has no "mass" term (a term with no gradients) in £0-1; the constraints on 
the couplings ensure that no such mass term is ever generated, and that the fluctuations of 
ip remain gapless. This is exactly what is expected as the Neel order parameter is non-zero 
in both the C and D phases. Finally, note that although we have only displayed Z„ for 
H = 0, the form of the coupling to a field in the z direction can be easily deduced from 
the gauge invariance arguments of Ref M . 
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The transition between phases C and D is controlled at tree- level by the sign of r4, and 
phase C is present for > 0. In this case, we can integrate out the massive field (f), order by 
order, and obtain an effective action for t/j. At tree-level, this effective action corresponds 
to putting (f) — in Ca- The ip fluctuations then represent the two real gapless spin wave 



modes of the Neel phase with dispersion 00 — ^K2/Kik. The loop corrections from the 
fluctuations will simply renormalize the values of and Ki. 

Finally, we turn to phase D which exists for < 0. In this case, acquires a mean value 



Ui 



(3.7) 



We have chosen (0) to be real, corresponding to choosing a definite direction for the com- 
ponents of hn orthogonal to n^. We now look at small fluctuations about this state by 
writing 







' — r4 



+ {(j)^ + i(f)y)/V2 



(3.8) 



where all flelds on the right hand side are real. The held (f)x represents changes in the 

magnitude of the condensate; its fluctuations arc therefore massive and can be neglected. 
We Fourier transform the remaining flelds, introduce the vector $ = [ipx, ipy, (py) and obtain 
from Ci the quadratic action: 



/ 



(27r)<^+i 2 
where the dynamical matrix M is 



-$i(-/c, -u)Mij{k, uj)^j{k, cu) 



(3.9) 



M{k,uj) 







Kiuj \ 

Kiu^ + K2k'^ 

/ 



(3.10) 



The zeros of the determinant of M give us the eigenfrequencies of the harmonic oscillations 
about the ground state in phase D. Analytically continuing to real frequencies, we obtain one 
linearly dispersing mode cu = \J K2/ Kik and one quadratically dispersing mode (at small k) 

uj — [K^K2) I {K\ -I- KxK'ik^)]^ . Although there were three real flelds to begin with, there 
are only two distinct modes, as a pair of them behave like canonically conjugate degrees of 
freedom. 



B. Model without V symmetry 

Now M 7^ 0, and a new class of terms will be permitted in the effective action. 
The analysis in this case is technically similar to the M = case; so we will be very 
brief. The effective action functional for the n,, and h,, now becomes: 
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Vn^Dh^ exp J df^x dr^Ci + £2) j 




2 V ■ 2 

2^8 



+ K,{Vn,){Vh,) + '-^hl+'^{hlf 



A = f + r,h,n, + ^{nlf + + ^-^{n.h.f (3.11) 

Note the presence of the additional couphngs K4 and r^; these couphngs were forbidden 
earher by the V symmetry. 

We can again look for static, spatially uniform, saddle points of C As before, this 
gives three different types of solutions, corresponding to the phases A, C, and D but none 
corresponding to the B phases. The low-lying excitation spectrum can again be found by 
defining a hybrid sigma model by imposing the constraints 

Uf^rif, = 1 ; Uf^hf, = c (3.12) 

where c is some constant. Clearly the only difference from the previous section is the non- 
zero (but constant) value of n^/i^. This is necessary because though the onset of the canted 
phase D is still signaled by the appearance of an expectation value of in a direction 
orthogonal to the mean direction of n^, the component of /i^ parallel to is non-zero in 
both the C (Neel) and D phases when M ^0 . 

We proceed as before and move well away from phase A, assume that has a well 
developed expectation value along the z direction, and replace the earlier parametrization 
by the following: 

( ip + ijj* ip — ip 
n = \ - 




The new feature is the first term in the second equation. From now on the analysis is similar 
to that in the previous section, so we merely state the results. At the Gaussian level, there 
now is a new term in the lagrangian C^i equal to K^{S/il)*'W(l) + VipVip*). Both the Neel 
and canted phases continue to have the same excitation spectrum as at M = 0, but the 
constants of proportionality in the low k dispersion relations are different from their M = 
values, as are the actual eigenvectors corresponding to the normal modes. 



IV. QUANTUM PHASE TRANSITIONS 

In this section we consider the critical behavior of some of the continuous quantum 
transitions among the phases in Fig |^. The transition between the paramagnet A and the 
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Neel phase C has aheady been discussed in some detail in Ref p2|. We will present below 



a discussion of the transitions from the Neel (C) and quantized ferromagnet (B) phases to 
the canted phase D. There are also a number of multicritical points in Fig|l|, but we will not 
discuss them here. As in previous sections, the discussion below will implicitly assume that 
d > 1, and we will obtain results for all such d. 



A. Neel (C) to Canted (D) Transition 

We begin with some general scaling ideas. This transition involves the onset of a mean 
value in a conserved charge, and this could possibly impose constraints on the critical ex- 
ponents. In an earlier paper 0, one of us had studied such constraints for the case of a 
transition from a paramagnet to a partially polarized ferromagnet. In the present case, the 
'paramagnet' is the Neel phase C, and the orientation of the Neel order parameter has As 
we show below, this is an important difference, which modifies the scaling relations. 

The order parameter describing this transition is the coarse-grained angular momentum 
density L±{x,t) in a direction perpendicular to the direction of Neel ordering which we 
choose to be the z axis. The response of the system to a field that couples to L± is described 
by the correlation function 

X^\k,uj) =1 J d'^xdt e^^{-i(k ■ X - u}t))i\L\x,t),L\Q,Q)])e{t) = 5^^xik,u), (4.1) 

with a, 6 = 1, 2, and where the angular brackets denote both thermal and quantum averaging. 
For simplicity, we will only discuss here the scaling at zero temperature. The behavior at 
finite temperature can be obtained from finite size scaling through standard arguments 



25,22]. On approaching the transition from the the Neel phase, we expect that x satisfies 



the scaling form 



X(fc,^) = -^^(fce,aiu;r) (4.2) 



where ^ is the diverging correlation length associated with the second order transition. 
g{x,y) is a universal function of it's variables and Oq and Oi are (possibly non-universal) 
constants. The susceptibility to a static, uniform external field that couples to L± is given 
by xo = linifc^o linioj^o x(^) ^) diverges with an exponent 7 = i^(2 — rj). It is possible to 
derive a scaling law relating z and rj, if we make the additional assumption that the stiffness, 
Ps, of the system to twists in the direction of the Neel order parameter has a finite non-zero 
value at the transition. While we have no proof that this is so in all cases, it is physically 
reasonable in view of the fact that the Neel order parameter is finite and non-vanishing 
across the transition. The explicit computations below will obey this assumption. 

Throughout the Neel phase at very low k,uj, x{k,uj) has the form predicted by hydro- 
dynamics: 

X{k,uj)= ^'''l ^ (4.3) 
with the spin wave velocity Vg given by 
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The divergence of xo at the transition thus imphes that the spin-wave velocity vanishes. 
In the vicinity of the critical point, we require that the scaling form ( [4.2|) reduce to the 
hydrodynamic expression ( [4.3| ) in the limit fc^ — >■ 0, aiuj^^ — > 0. For this to happen, the 
function g{x,y) must satisfy 

9ix,y)-*— as x,y^O (4.5) 

— 

This gives the scaling of xo and v"^ as xo ~ ■C^"^, v'^ ~ ^2(1-2) _ (j^he xo scaling is equivalent 
to the relation 7 = z/(2 — r^).) From the hydrodynamic expression ( [4.4| ), we then get 

z = 2 - I (4.6) 

Below we will verify this scaling relation by an explicit renormalization group (RG) calcu- 
lation. This scaling relation replaces the relation z = d + 2 — t] proposed earlier [0 for the 
case of onset of ferromagnetic order from rotationally invariant paramagnetic phases; as we 
will see in more detail below, the latter relation has been modified here because 1/ps (the 
dimensionful inverse of the stiffness of the Neel order) behaves like a "dangerously irrelevant 
coupling" . 



1. Renormalization group analysis; V symmetry present 



We will examine here the behavior of the action (|3.6|) for systems with V symmetry 
under RG transformations. As noted earlier, this action undergoes a phase transition as 
a function of the tuning parameter from a Neel phase C (r4 > in mean-field theory) 
with {(f)) = to a canted phase D (r4 < in mean-field theory) with (0) 7^ 0. The field ip 
represents fluctuations of the Neel order about its mean value and we always have (ip) = 0. 

Notice that the ip field plays a role analogous to the spin- wave modes of the usual 0(3) 
non-linear sigma model |^ : their correlations are always gapless, and no mass-term for them 



is ever generated. It therefore seems natural at first to construct an RG under which the field 
ijj remains dimensionless and its effective stiffness K2 fiows. So for the RG transformation 

x' = x/s t' = tIs\ (4.7) 

where z is an unknown exponent, we have ip'{x',T') = ip{x,T). At tree level invariance of 
the action requires 

K'^ = K2s'^+~'~^ (4.8) 

In contrast, the field is a soft-spin order parameter, and should behave under RG like 
a conventional Landau-Ginzburg field. In this case, we fix K3 = constant which requires 

0' = (4.9) 
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at tree level. Now demanding that the Ki{(f)*dr'ip — (pdrip*) be invariant we obtain, again at 
tree level 



z = d + 2 (4.10) 

This appears to be the tree-level version of the identity z = d + 2 — rj discussed in Ref 
for the ferromagnetic onset transition, and it is tempting therefore to identify z as the 
dynamic critical exponent. Notice, however that for this value of z, K2 flows to infinity for 
all (i > 0, and there is thus the possibility that the dimensionful constant could appear 
in determining an important low-frequency scale which is not characterized by the exponent 
z. This is indeed what happens, as we can easily see by studying the Gaussian fluctuations 
at the critical point r4 = 0: the fluctuations involve modes at a frequency uj ~ J K2K3/ Klk"^ 



with momenta k. This dispersion is consistent with the scaling ( |4.(j| ) of K2 and the value 
( [4.10| ) for z. However the present approach is quite awkward, and it is clear that a scheme 
with dynamic exponent z = 2 would be preferable. 

Such a scheme is motivated by the realization that ^ 00 in the above RG (and so 
I/K2 is dangerously irrelevant). As a large K2 suppresses fluctuations of ip, we can safely 
integrate out the ip field completely from ( |3.6| ), with the expectation that only the leading 
terms in I/K2 will be important. Such a computation leads to the following effective action 
for the field 

S = J + + \<Pik,-^)n^Jd'xdr\<p{x,r)\\ (4.11) 

We have rescaled the field (p/y/K^, and introduced the couplings 60 = Kl/{K2K^)^ 
tq = r4^/K^ and u = Ui/K^. Notice now that there is a long-range term associated with the 
uj'^/k'^ non-analyticity in the propagator. Standard methods to analyze field theories with 
long-range couplings and anisotropic scaling in spacetime are available and they can 



be applied to the action S in (^4.11|) in a straightforward manner. The non-linearities in S 



are controlled by the dimensionless coupling constant 

go = ub,'/'sy-^ (4.12) 

where Sd = 27r'^/V(r(rf/2)(27r)'^) is a phase space factor and /i is a momentum scale. The 
dependence on n indicates that the Gaussian fixed point controls the infrared behavior for 
d > 2: this fixed-point has the exponents z = 2, rj = 0, and u = 1/2. 

For d < 2, a. RG analysis using an expansion in powers of e = 2 — ci is necessary. As usual, 
we define the renormalized field = A/Z0ij, and the renormalized couplings go = {Z/i/Z'^)g 
and 60 = {Zb/Z)b. Note that the coupling g here should not be confused with the g in the 
quantum rotor Hamiltonian (|1.3| ). The renormalization constants Z, Z/^^ Z}y are determined 
in a minimal subtraction scheme and we obtained (some details appear in Appendix y): 

Z = l-^ + ... 
72e 

5(7 

Z,^l + ^ + ... 

Zb = l (4.13) 
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The result for Z\, is expected to be exact to all orders in g — there are no uJ^jk^ terms 
generated in the two-point function in the loop expansion. The above results required the 
evaluation of certain one and two-loop Feynman graphs — these are most easily evaluated 
by writing each propagator in spatial and time co-ordinates and then doing the necessary 
integrals (Appendix |C|). We can now obtain the /3 function for g: 



B 



-e9 + — 



(4.14) 



where the derivative is taken for a fixed bare theory. The f3 function has a fixed point at 
g = g* = 4e/5, which is the infrared stable fixed point for d < 2. Scaling exponents of the 
critical point can now be determined: 



^ 1 7 

a/i 



1 d Zb 



4e^ 
225 



B,a=g* 



2 - 



2e2 
225 



(4.15) 



Note that the scaling relation ( [4.6|) is obeyed, and is a consequence of = 1. The value 
of r controls the deviation of the system from criticality, and the renormalization of cjP' 
insertions will determine the value of the critical exponent v. We define its renormalization 
by ro = {Z2/Z)r, and obtain in a similar manner 



1 + ^ 
2e 



(4.16) 



which leads to 



1 . d Z2 
- = 2 - 11— In — 
u oil Z 



B,g=g* 



5 



(4.17) 



2. Renormalization group analysis; V symmetry absent 



As we noted earlier in Section [llllj| , in the absence of V symmetry, the most important 
modification of the action (|3^ ) is that a term K^{{y ifj*) .{V (p) + (VV')-(V0*)) is now per- 
mitted. As in Section |1V A 1| we may now integrate out from such an action, which leads 



to the following modified version of (|4.11 



S 



d'^k doj 
(27rp27r 



i\Qh]!^uj + ro\\(t){k,uj)\' + ^ I d''xdT\(f){x,T)\\ (4.18) 



u 



k^ 



The new feature is the iu term in the propagator. Associated with this term is the dimen- 
sionless coupling Aq; power-counting shows that Aq marginal at tree level. Indeed Aq remains 
marginal at one- loop, and a two- loop calculation is necessary to decide if the Aq = fixed 
point of Section [IV A 1| is stable. 
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We define the renormalization of A by XoBq^'^ = [Zx/ Z)\h^/'^ . Tlie new values of tlie 
renormalization constants can now be computed to be (Appendix y): 



1 - 

1 + 
1 

1 + 
1 - 



A2 + 9 



2e(A2 + 4)(2A2 + 9)2 
g A2 + 20 



+ . . . 



2e (A2 + 4)3/2 

4^ 1 

e (A2 + 4)3/2 

9' 1 



+ 



+ . . . 



2e(A2 + 4)(2A2 + 9) 



+ . . . 



(4.19) 



Tliese results require an expansion in powers of g, but the dependence on A is exact at each 
order. The j3 functions of g and A can now be determined 



R ^ 



~eg + 



B 



^2 A2 + 20 
2 (A2 + 4)3/2 
g^ A(5A2 + 27) 



(4.20) 



B 



2 (A2 + 4)(2A2 + 9)2 

One of the fixed points of ( [4.20| ) is the model with V symmetry of Section |IV A 1 

A* = 0, g* = 4e/5. (4.21) 
However, this fixed point is unstable in the infrared; the stable fixed point is instead 



A* 



00, 



(4.22) 



for A ^ 00, the loop corrections become functions of g/X, which is therefore the appropriate 
measure of the strength of the non-linearities. Near this stable fixed point, the frequency 
dependence of the propagator is dominated by the iu term and the uj'^/k'^ term can be 
neglected. The exponents at this fixed point can be determined as in ( [4.15|) and ( [4.17 ) from 
( [4.19| ) and we find the same exponents as those of the d > 2 Gaussian fixed point: 



77 = 



V 



1/2 



(4.23) 



Although these exponents are the same for d > 2 and for d < 2, it is important to note 
that the d < 2 theory is not Gaussian — there is a finite fixed point interaction (g/X)* = 2e 
which will make the scaling functions and finite temperature properties different from those 
of the Gaussian theory. The structure of this non-Gaussian fixed point has been explored 
earlier pSj^ in in a different context and more details are provided in Appendix 0). Also 
note that the flow of A to this fixed point is slow, and will lead to logarithmic corrections to 
scahng. 
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B. Quantized Ferromagnet (B) to Canted (D) Transition 



The ferromagnetic order parameter (L^) is non-zero in both the B and D phase, while 
the component of orthogonal to (L^) plays the role of the order parameter. The roles of 
and are therefore reversed from the considerations in Sections |T| and [IV A| , and the 



analysis here will, in some sense, be dual to the earlier analyses. In the following we will 
assume that M 7^ and that V is not a good symmetry; the special properties of the M = 
case will be noted as asides. 

To begin, we need a rotationally-invariant field theory which describes the dynamics in 
one of the quantized B phases and across its phase boundary to the D phase. Let us work in 
the Bfp phase i.e. the phase in which the quantized moment is £0 7^ per rotor. If we now 
proceed with a derivation of the path integral as in Section |II|, we find that the single site 



functional integral (|3.1|) has to be evaluated in a regime in which ( p.2[ ) has a minimum at 
£ = Iq-. this is quite difficult to do. We circumvent this difficulty by the following stratagem. 
We consider a somewhat different model, which nevertheless has a quantized ferromagnetic 
phase (with a moment £0 per site) and a canted phase separated by a second order transition; 
this model has an enlarged Hilbert space, but the additional states have a finite energy, and 
we therefore expect that its transition is in the same universality class as in the original 
rotor model. The new model is an ordinary Heisenberg ferromagnet coupled to a quantum 
rotor model through short-ranged interactions. The Hamiltonian Hnew = Hfjeis + Hrot + Hint 
where 

Hint = -(E Ki-jSi^U^ + E Mi_jSi^n,^) (4.24) 

Here Si^ is an ordinary Heisenberg spin operator acting on a spin £0 representation, and it 
commutes with the and L^. The Ki_j and Mi_j are some short-ranged interactions and 
A is a positive constant. (There is an overlap in the symbols for the couplings above with 
those used before; this is simply to prevent a proliferation of symbols, and it is understood 
that there is no relationship between the values in the two cases.) The total conserved 
angular momentum is J2i + Si^. The Hilbert space on each site consists of the spin £q 
states of the Heisenberg spin and the £ = 0, 1, 2, 3, . . . states of each rotor. Ignoring inter-site 
couplings, the lowest energy states on each site are a multiplet of angular momentum io, 
similar to the B^^ phase of the original model. 

It is now easy to see that Hmw exhibits a quantized ferromagnetic phase (with a moment 
io per site) and a canted phase for appropriate choices of parameters. First imagine that 
all the Ki_jS and Mi_j^s are zero and that g ^ J. Then naive perturbation theory in J/g 
tells us that the ground state is a quantized ferromagnet with total angular momentum io 
per site. Now imagine introducing non-zero values of the K and M couplings. The rotor 
variables then see an effective "magnetic field" that couples to the Lj^, and an effective field 
that couples to the hi^. The latter has the effect of introducing an expectation value for the 
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component of parallel to (S*^). The "magnetic field" does not introduce an expectation 
value of if it is weak. As the strength of the "magnetic field" is increased by varying the 
Ki-j^s, it induces a transition into a phase with a non-zero, non-integer expectation value of 
L^, and a non-zero expectation value for the component of perpendicular to (S^). This 
is the transition to the canted phase. 

The derivation of the path integral of this model is easy; because the Hilbert space 
has now been expanded, we have also have to introduce coherent states of spin io at each 
time slice in the Trotter product. The final functional integral will be expressed in terms 
of the field of the quantum rotor, and a unit vector field (s^ = 1) representing the 

orientation of the spin S^. Finally, as in Section |T|, we also impose the constraint n^s^ = c 
to project out fluctuations which have a local energy cost (this cost is due to the effective 
field, noted above, that acts between the and ra^). Our final form of the modified sigma 
model appropriate for the transition from a quantized ferromagnet Bi^ to the canted phase 
D is 

Z = J Vn^VSf,5 (s^ - l) (5 (n^s^ - c) exp ^- j <fx drC^ 



Os I On 

+ ^iyn.f + ^(V.,)^ + K,{Vn,){Vs,) + + "^{nlf (4.25) 




Again there is an overlap in the symbols for the couplings above with those used in Sees |T1 



and [IV B| , and it is understood that there is no relationship between the values in the two 
cases. There is no length constraint on the components of orthogonal to s^, which act 
as a soft-spin order parameter for the transition from phase B to phase D. The pSy term is 
the effective internal "magnetic field" (noted above) that the Heisenberg spins impose on 
the n^, with p a coupling constant that measures its strength. Note also that 

= c = if P is a symmetry. (4.26) 

The coupling Mq = loa~'^, where a'^ is the volume per rotor, and, as before, is the external 
magnetic field. The first term in C is the Berry phase associated with quantum fluctuations 
of the ferromagnetic moment: ^^(s) is the vector potential of a unit Dirac monopole at the 
origin of s space, and is defined up to gauge transformations by ef^u\dA\/dsu = s^. 

The onset of phase D is signalled by a mean value in orthogonal to the average 
direction of s^. Using a reasoning similar to that above ( p.5|) , we parametrize the fields of 
the modified sigma model as 

2 (2-|0^v^i^(2-|0|Y^l-l0^) 

_ f ^+r i^-r (^-0+r^)(2-i0p)v^ \ 
1 V2 ' 72^ ' 72(1 -m ^"-''^ 
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We have chosen a non-standard parametrization of the unit vector field in terms of 
a complex scalar 0: this choice is a functional integral version of the Holstein-Primakoff 
decomposition of spin operators, and ensures that the Berry phase term iMoA^(s)(9s^/(9r) 
is exactly equal to Mo0*((90/9r). The parametrization of then follows as in (|3.5| ) and 
ensures that n^s^ = c. The two complex fields ip, 0, correspond to the 4 degrees of freedom 
of the hybrid sigma model. The field (f) represents oscillations of the s about its mean value: 
by definition, we will have (0) = in both the B and D phases. The field if) measures 
the amplitude of orthogonal to the instantaneous value of s^: we have {tp) = in the 
quantized ferromagnet B, while {ip) 7^ in the canted phase D. We now insert ( [4.27 ) into 
(|]2|) and obtain 

OT OT \ OT OT J 

+ (-Mo + + K^\iP\^ + K,{rct> + W)) 

+ K2I VV^P + + K^{VrS/(l) + VV^V0*) 

+ ^2(1^1^ + \{r^ + ^^*?) + y IV^r (4.28) 

where = 2pKi, Kq = cpKi, Kj = K3 + 2cKi + c^i^a, = + CK2, = - 
Kip^ + MiC^/2. Stability of the action requires that K2 > 0, K-i > 0, and > K^] these 

conditions will be implicitly assumed below, and are the analog of the requirement in the 
mean-field theory of Section y that M not be too large. Without loss of generality, we can 
assume that Mq > 0; however the signs of and MqK^ — Kg can be arbitrary, and these 
will play an important role in our analysis below. Systems with M = 0, and the symmetry 
V, will have Kq = Kg = 0. We have explicitly displayed only terms up to quartic order in 
the fields. Anticipating the RG analysis below, we have also dropped irrelevant terms with 
additional time derivatives. The field if^ = (0, 0, H^) has been assumed to point along the 
z direction. 

Notice that the coupling to could have been deduced by the gauge invariance prin- 
ciple which demands the replacement d/dT d/dT + [d/dT d/dT — Hz) when 
acting on the 0, (0*, ip*) fields. The magnetization density, (Ai) is given by 

(M) = -§-, (4.29) 

where is the free energy density. The gauge invariance of the action ensures that at T = 0, 
this magnetization density is exactly Mq, as long as (-0) = i.e. we are in phase B. 

We can easily do a small fiuctuation, normal mode, analysis in both phases of ( [4.28|) . In 
the quantized ferromagnetic phase B (r2 > 0) we have the usual ferromagnetic spin waves 
with uj = Kik'^/Mo at small k and also a massive mode with u ~ r2/{K^ — Kq/Mq) (we 
have assumed here and in the remainder of the paragraph that Hz = 0). The massive mode 
becomes gapless at r2 = 0, signaling the onset of the canted phase D. In this canted phase. 



we have (-0) = y—r2/ui. Analysis of fluctuations about this mean value for -0 gives two 
eigenfrequencies which vanish as A; — > 0; the flrst has a linear dispersion 
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while the second disperses quadratically 



W 



(4.31) 



where W = (Mq/^ - K^^if + 2KqMq{^/KJCj - Kg) > 0. These results for phase D 
are in agreement with those obtained earlier using a 'dual' approach with the action ( p.6|) 
in Sections |111 A| and [111 Ij . 

for the transition from C to D 



It is also useful to compare and contrast the action 
with the action ( [4.28| ) above for the B to D transition. Notice that the roles of the fields 
are exactly reversed-in the former was the order parameter for the transition while ip 
described spectator modes which were required by symmetry to be "massless" ; in the latter 
ip is the order parameter while is "massless" spin-wave mode. The two models however 
differ significantly in the nature of the time derivative terms in the action. In ( [4.28|) there 
are terms linear in time derivatives for both the ip and fields, while there are no such 
terms in ( p.6|) (the time derivative term involving an off-diagonal 0, ip coupling is present 
in both cases, though). As a result the spectator spin-wave mode in ( ^4.28| ) has u ~ k"^, 
while the spectator spin-wave ip mode in (|3.6| ) has ~ A;. These differences have significant 
consequences for the RG analysis of ( [4. 28] ), which will be discussed now. 



1. Renormalization group analysis 
The logic of the analysis is very similar to the 'dual' analysis in Sections [IV A 1| and [IV A2. 



We expect the stiffness of the spectator modes flows to infinity, and hence it is valid to 
simply integrate out the fluctuations. This gives us the following effective action for the ip 
field (which, recall, measures the value of in the direction orthogonal to the ferromagnetic 
moment): 

[ d'^k du ( .^^ ^^,2 i-iKeUj + Ksk'^f \,,n m2 
^ = y (2^ 2^ (-^^^^ + - -^M,. + KJ + j I 



-fyl d''xdT\^{x,T)\\ (4.32) 



We have written the action for Hz = 0, and the Hz dependence can be deduced from the 
mapping —iu —iuo + Hz- Simple power-counting arguments do not permit any further 
simplifications to be made to this action. Despite the apparent formidable complexity of its 
form, the RG properties of ( 4.32[ ) are quite simple, as we shall now describe. 



We will consider the cases with and without V symmetry separately. 



a. V symmetry present 

Now we have Kq = Kg = 0. The non-analytic terms in the propagator of ([4.32[) disappear, 
and the remaining terms are identical to those in the action of a dilute Bose gas with a 
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repulsive interaction ui. The order parameter, and the spectator modes, 0, are essentially 
independent, with all couplings between them contributing only "irrelevant" corrections to 
the leading critical behavior. 

The critical theory of the quantum transition in the dilute Bose gas has been studied ear- 
lier p5| , p9| (see Appendix 0), and also appeared as the fixed point theory in Section [IV A2 



A measure of the strength of the loop corrections is the dimensionless coupling constant qq 
(the analog here of Eqn ( [4.12| )) 

go = uiK:,^\K^\-^Sd^''~^ (4.33) 

(where, as before, is a renormalization momentum scale, Sa is a phase space factor, and 
the symbol g is not to be confused with the coupling g in ( |1.3| )), which determines the 
scattering amplitude of two pre-existing excitations. This amplitude undergoes a one-loop 
renormalization due to diagrams associated with repeated scattering of the two excitations, 
a process which leads to the (5 function 

^ = -^9+^^ (4.34) 

where e = 2 — d and we have used the minimal subtraction method to define the renormalized 
g (this (3 function is in fact exact to all orders in g — see Appendix 0). The infrared stable 
fixed point is g* = for d > 2, and g* = 2e for d < 2. The exponents take the same values 
as in Section IIV A 2 



z = 2 ri = v= 1/2. (4.35) 
for all values of d. It is important to note, however, that the present subsection and Sec- 



tion |1V A2| have very different interpretations of the order parameter. In Section [IV A 2| the 
exponent t] is associated with the field which is proportional to the ferromagnetic moment. 
Here the r] refers to the the field ip, and the ferromagnetic order parameter scales as as 
we discuss below. 

The magnetization density order parameter A^(x,r), is defined by generalizing (|4.29[ ) 
a space-time dependent external magnetic field Hz{x,t). In the quantized ferromagnetic 
phase B, clearly (A^) = Mq. Therefore, the deviation Mq — M. can serve as another order 
parameter for the B to D transition. Indeed, this order parameter scales as I'i/'p, and we can 
define an associated "anomalous" dimension rjM- A simple calculation shows that rjM = 
and so the identity 

z = d + 2-riM (4.36) 



is always satisfied. As in Ref [^, we can also compute the behavior of the mean value of 



this conserved order parameter in the canted phase D. For d < 2, this behavior is universal 
and is given by 

{M) = Mo-sgn{K,)e{-r2)cJg^] , (4.37) 



2 
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where Cd is a universal number; we compute Cd in an expansion in e in Appendix 0. Note that 
the magnetization density can increase or decrease into the ferromagnetic phase, depending 
upon the sign of K^; this behavior is consistent with that of the mean field theory in 
Section ||. Note also that ([4.37|) is independent of the magnitude of K^. The magnitude 
of K5 however does determine the width of the critical region within which ( [4.37| ) is valid. 



For IK^I smaller than ^/r2, it becomes necessary to include the leading irrelevant frequency 
dependence in the propagator of ( [4.32| ) - a term. The point i^'s = is a special 

multicritical point at which this new term is important at all values of this multicritical 
point has z = 1 pH|, and has {A4) continue to equal Mq (to leading order) in phase D. 



b. V symmetry absent 

The theory without V symmetry is somewhat more involved, and it is convenient to work 
with a more compact notation. We rescale ip{k,u) {Mo/^/WK^)■^|J{k, u), uj ojK-j, and 
manipulate ( 4.32 ) into the form 



+ ^ J d'^xdT\^{x,T)\\ (4.38) 

where ro = M^r2/W, u = uiKjM^/W^, Aq = Kj{MqK^ - Kl)/W, and 60 = (^6^7 - 
K^MqY /W K^. Notice that we have introduced two dimensionless couplings, Aq and 
while Mq has the dimensions of time/ (length)^. The stability conditions discussed below 
i WM ) imply that < 60 < 1- 

We now discuss the renormalization of the theory ( [4. 381 ) as a perturbation theory in the 
dimensionless coupling constant 

r?o = «5,/-V(Mo|Ao|), (4.39) 

which is the analog of ( [4.33| ). A key property of this perturbation theory, similar to that 
noted in Section IV A 1| , is that no terms sensitive to the ratio w/fc^, as u; — > 0, — >■ 0, are 



ever generated. As a result the term multiplying 60 in (|4.38D does not undergo any direct 



renormalization. With this in mind, we define the renormalized field hy = ^fZipn^ and the 
renormalized couplings by Aq = {ZxlZ)\, bo = {1/Z)b, ro = {Z2/Z)r and go = {Z4^/Z'^)g. It 
remains to compute the renormalization constants Z, Zx, Z2, and Z4. 

In the perturbation theory in g, the sign of A (which is also the sign of MoK^ ) plays 
a crucial role. The physical interpretation of this sign is the following: when A > 0, the 
spin wave quanta, 0, and the order parameter quanta, ip, carry the same magnetic moment; 
however for A < they carry opposite magnetic moments. This was already clear in the 
Kq = result (|4.37| ) , where the condensation of lead to a decrease or increase in the mean 
magnetization density depending upon the sign of K^,. For A < 0, it is then possible to have 
a low-lying excitation of atp and a quantum, which has the same spin as the ground state. 
For the system without V symmetry, such an excitation will mix with the ground state, and 
leads to some interesting structure in the renormalization group. In the actual computation, 
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the sign of A determines the location of the poles of the propagator of (|4.38|) in the complex 
frequency plane. A simple calculation shows that the propagator has two poles, and both 
poles are in the same half-plane for A > 0, while for A < the poles lie on opposite sides of 
the real frequency axis. We will consider these two cases separately below: 
A > 

For the case where the poles are in the same half-plane, we can close the contour of frequency 
integrals in the other half-plane, and as a result many Feynman diagrams are exactly zero. 
In particular, all graphs in the expansion of the self energy vanish. As a result, we have 

Z = Zx = Z2 = \ (4.40) 

There is still a non-trivial contribution to Z4 (from the 'particle-particle' graph): 

9 (A + 1-6) 



^'"^+2e(l-6)(l + A) 
where, as before e = 2 — (i. These renormalization constants lead to the /? functions 



(4.41) 








(l + A-6) 
2 (1-6)(1 + A) 



(4.42) 



The vanishing of and /5a is expected to hold to all orders in (7; as a result, h and A 
are dimensionless constants which can modify the scaling properties. The exponents are 
however independent of the values of h and A, and are identical to those of the system with 
V symmetry, discussed above in Section p.V B l| a. The upper critical dimension is 2, and for 
d <2 the magnetization density obeys an equation similar to ([4.37|) 



l^M) = Mo - 0(-ro)C,|ro|"/^ 



(4.43) 



The 'universal' number Cd now does depend upon b and A: the leading term in Cd can be 
computed by using the gauge-invariance argument to deduce the modification of ( [4. 381 ) in 
a field, using ( |4.29| ) to determine the expression for A4 [^], determining the fixed point of 
j3g, and then using the method of Ref ||29l 



-this gives 



Cd 



6)(1 + A) 




(4.44) 



This result should be compared with the first term in ( pi2 ), which is the value of Cd for the 
Bose gas. Terms higher order in e are also expected to be modified by b and A, but have 
not been computed. The results (OD, ( ^441) are not useful for too small a value of A: the 
reasons are the same as those discussed below ( [4.37]) but with A playing the role of K^. 
A < 

The computations for A < are considerably more involved and are summarized in Ap- 
pendix |^. From (p5|), (p6|) and ( p^ ) we can deduce the renormalization constants 
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Z=l-|Ji?i(6,A) 

^A = l-^i?2(&,A) 
Z,= l + j-R,{b,X) 

Z2 = l + ^R4{b,X), (4.45) 



where -Ri, R2, R3, Ra are functions of h and A defined in ( ^Tj^Ol) . These now lead to the j3 
functions 

f3, = g%Ri{h,\) 

f3x = 9\XRi{b,X)-R2{b,X)) 

Pg = -eg + g^Rsib,X)/2. (4.46) 

It is not difficult to deduce the consequences of these flows by some numerical analysis aided 
by asymptotic analytical computations. 

For c? > 2 (e < 0) we have 6(/i) b*, A(/i) A*, g{^) in the infrared (/i 0). 
Here < 6* < 1 and A* < 0, but the values of b*, X* are otherwise arbitrary and determined 
by the initial conditions of the flow; they acquire only a finite renormalization from their 
initial vales. The fiow to this final state has a power-law dependence on /i for e < 0, while 
it behaves likes ~ l/log(l//i) for e = 0. The properties of the final state have only minor 
differences from A > case discussed above, and we will not elaborate on them. 

The possible behaviors are somewhat richer for e > 0. The flow in the infrared is either 
to a flxed line or a flxed point. 

(z) The flxed line is 2/3 < 6* < 1, A* = 0, (?* = 2e; the flow of g to g* is a power-law in fi, 
while that of b and A is logarithmic: 



b*) 



L*\2 



^^^^ 2e26*(36*-2)log(l//i) ' ^^^^ ^* 4e2(36* - 2)2 log(l//i) ' ^^'^'^^ 

(ii) The flxed point is b* = 0, A* = —00, g* = 2e; again the flow of g to g* is a power-law in 
II, while that of b and A is logarithmic: 

A(^) ~ -(log(l//i))^/^ , %) ~ (log(l//i))-V3. (4.48) 



The exponents at both the flxed line and the flxed point take the same values as in ( 4.35|) . 



However the flows above will lead to logarithmic corrections which can be computed by 
standard methods. The prefactors of the these corrections, and also some amplitude ratios, 
will vary continuously with the value of b* along the flxed line. 



V. ROTORS WITH A NON-ZERO MIMIMUM ANGULAR MOMENTUM 

We now consider an extension of the rotor model H in which each rotor has a 'magnetic 
monopole' at the origin of n space. Our motivations for doing this are: (i) It allows for 
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quantized ferromagnetic states in which the moment is a half-integer times the number of 
rotor sites. Such phases will occur in cases in which each rotor corresponds to an odd number 
of Heisenberg spins in an underlying spin model, {ii) In d = 1 (to be discussed in Sec |VID 



such rotor models will lead, by the construction of Ref fSTJ, to Neel phases described by a 
sigma model with a topological theta term. 

Following the method of Wu and Yang , we generalize the form (|1.2| ) of to 

Lif, = -e^yxriiy {^-^ + qSiAxini^)^ - qSirii^ (5.1) 

where q is chosen to have one of the values 1/2, 1,3/2,2,5/2, . . . and Afj_{n) is the vector 
potential of a Dirac monopole at the origin of n space which satisfies e^i,xdAx/dnu = (we 
considered this same function in Sec [IV B| for different reasons) . We will work exclusively 
on bipartite lattices and choose Ei = 1 {si = —1) on the first (second) sublattice: we will 



comment below on the reason for this choice. It can now be verified that (|5.1|) continues to 



satisfy the commutation relations (1.1) for all q. However the Hilbert space on each site is 



restricted to states \i, m) satisfying m) = i{i + l)|g, £, m) with £ = g, g + 1, g + 2, . . . 

and m = —i, + Notice that there is a minimum value, q, to the allowed 

angular momentum. 

Another consequence of a non-zero q is that it is no longer possible to have P as a 



symmetry of a rotor Hamiltonian. There is a part of in (5J.) which is proportional to 
n^, and this constrains L^, to have the same signature under discrete transformations: 
this rules out P as a symmetry, even in models with M = 0. 

Our insertion of the Ei in ( |5.1D is also related to the absence of the V symmetry for g > 0. 
Because and have the same signature under all allowed symmetries for g > 0, their 
expectation values turn out to be proportional to each other on a given site, with the sign of 
the proportionality constant determined by Ei (we will see this explicitly below). In a g > 
model with = 1 on every site, the spatial average of (n^) is proportional to that of (L^) 
and therefore to the net ferromagnetic moment; as a result, such models turn out to have 
only quantized ferromagnetic phases. Staggering of is a way of inducing a Neel-hke order 
parameter with no net moment; then a spatially uniform value of (n^) represents a staggered 
mean value of the angular momentum (L^), as expected from a Neel order parameter. The 



staggering of Ei was also present in the d = 1 analysis of Ref |]3T|, and was crucial there in 
generating the topological theta term. 

We present in Fig ^ the mean field phase diagram of the Hamiltonian H (Eqn ( |1.3| )) 
on a bipartite lattice with M = 0, given by (|5.1|) , and g = 1/2. Because of absence of 
the V symmetry, a non-zero M will not make a qualitative difference. Other values of the 
g are also expected to have a similar phase diagram. The caption describes the non-zero 
components of (L^) and (n^) in the various phases. The expectation values (n^) have exactly 
the same form as the (L^) in spin space, but the opposite signature in sublattice space {e.g. 
a staggered configuration of (L^) implies a uniform configuration of (n^) and vice versa). 

The phases in Fig || are closely related to those in Fig 0. As before we have the quantized 
ferromagnetic phases B (albeit, now with half-integral moments), the Neel phase C, and the 
canted phase D. The main difference is in the absence of a quantum paramagnetic phase 
A: this is clearly due to the minimum allowed value g = 1/2 in the single rotor angular 
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momentum. 

It is now possible to undertake an analysis of the low energy properties of the phases, 
and of the critical properties of the phase transitions, much like that carried out in Sees ^ 
and 1^ for the q = case: such an analysis shows essentially no differences between q = 
and g > 1, at least for d > 1. The excitations of the phases B, C, and D in Fig ^ are 
the same as those of the corresponding phases in Fig |l], as are the universality classes of 
the continuous phase transitions between C and D, and between B and D. One small, but 
important, point has to be kept in mind in this regard. The g > models are not invariant 
under the symmetry V even when M = 0, and so the restrictions that the V symmetry 
implies for the g = analysis of Sees |T| and |V| must not be imposed now. 

In d = 1, there are significant differences between the q = and q = 1/2 cases, and these 
will be discussed in the next section. 



VI. QUANTUM ROTORS IN ONE DIMENSION 

The general topology of the mean-field phase diagrams Figs |l| and ^ is expected to 
hold for all d > 1. As we have seen in Sec for d < 2 fluctuations do modify the 
critical properties of the continuous transitions, but these modifications are computable in 
a systematic expansion in 2 — d. In c? = 1, fluctuations modify not only the transitions, but 
also the stability of the phases: as a result, we expect significant changes in the topology of 
the phase diagram itself. Further, we also expect a sensitive dependence to the value of the 
monopole charge q. In the following, we present a mixture of results, educated surmises and 
speculation on the nature of the phases and phase transitions in d = 1 for different values 
of q. 



A. q = 

The expected phase diagram is shown in Fig ^. We discuss some of the important fea- 
tures in turn: 

(i) There is no phase with Neel order (the analog of phase C for d > 1); it has been pre- 
empted mostly by the quantum paramagnet A. Fluctuations in the incipient Neel phase 
would be described by a 0(3) non-linear sigma model in 1 -|- 1 dimensions (without a topo- 
logical term), which is known not to have a phase with long-range order. 
(m) The quantized ferromagnetic phases, B, are stable even in d = 1. They have the usual 

~ A;^ spectrum of spin-wave excitations. 
{in) The canted phase D has been replaced by a novel new phase E - the partially polarized 
ferromagnet. The phase E has true long-range order in the ferromagnetic order parameter 
(Lz); however the value of (L^) is not quantized and varies continuously. The canted phase 
D,in d > 1, also had long-range order in n in the x — y plane; in contrast, in phase E {d = 1) 
this long- range order has been replaced by a quasi- long-range XY order i. e. correlations of 
n decay algebraically in the x — y plane. The only true broken symmetry is that associated 
with (Lz), and the system is invariant with respect to rotations about the z axis (compare 
with phase D in which the 0(3) symmetry was completely broken and there was no invariant 
axis). Note that though there is no LRO in the x — y plane in phase E, the Goldstone mode 
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of phase D, associated with it's x — y LRO survives i.e. phase E has not only a gapless mode 
ui ^ k"^ due to ferromagnetic long-range order, but also a gapless mode ui ^ k due to quasi 
LRO in the h field. All of these results follow from a straightforward analysis of the actions 
i WM ) and (|;28D ind=l. 

{iv) A few remarks about the universality class of the phase transition between phases B and 
E. The theory (|4.28|) should continue to apply. The expansion in e = 2 — (Section [IV B 1|) 
should be valid all the way down to = 1 although it may not be quantitatively accu- 
rate. The theory with V symmetry has critical properties identical to a dilute Bose gas, and 
was discussed in Ref and solved there by a fermionization trick. The theory without V 
symmetry, has non-analytic terms in the action, and is probably not amenable to a simple 
solution by fermionization. 

(v) The phase transition between phase A and phase E is always expected to be first-order. 
The quantum paramagnetic phase A has no gapless excitations and a vanishing spin sus- 
ceptibility. It is then difficult to conceive of a mechanism which could lead to a continuous 
condensation of the field /i^ in an action like (|3.3|) and (|3.6|) (or (|3.11|) for M 7^ 0): integrat- 
ing out the A^^ fiuctuations does not yield a negative contribution to the "mass", ri, of the 
field, as the spin susceptibility is zero. 



B. Integers q > 

The phase diagram is shown in Fig ^. All of the phases are identical to those discussed 
above for g = 0. The only difference is in the behavior of the first-order line surrounding 
phase A: it bends down towards the origin of the J-K plane, as the minimum possible value 
of the single site angular momentum always forces in quantized ferromagnetic phases for J 
small enough. 



C. Half- integers q > 

The phase diagram is shown in Fig |^. The topology is now similar to the positive integer 
case, as are the phases B and E, and the transitions between them. The primary difference 
is that the quantum paramagnet A has been replaced by a critical phase F, which has no 
broken symmetries and power-law decay of all observables. Fluctuations in phase F are 
described by the 1 + 1 dimensional 0(3) non- linear sigma model with a topological term 
with coefficient 9 = -k: this mapping follows from the analysis of Ref |^ . 



Finally, we make a few remarks on the transition between the critical phase F and 
the partially polarized ferromagnet E. As F has gapless excitations, this transition can be 
continuous. A (strongly-coupled) field theory for this transition is given by the action £1 
in (|3.11|) , supplemented by the constraints ( p.l2| ) and a topological term at 6* = tt in the 



unit-vector A^"^ field. An alternative, Hamiltonian point-of-view on the same transition is the 
following. It is known that the critical 1 + 1 dimensional 0(3) non-linear sigma model 
is equivalent to the SU{2) Wess-Zumino-Witten model [^] at level k = 1. The Hamiltonian 
of this model is given by 



Hwzw = (yiJl + jl) (6.1) 
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— * 

where Jl^r are the SU{2) currents obeying a Kac-Moody algebra. We now want to induce 
a ferromagnetic moment into the ground state of this theory. The action does this by 



couphng in a fluctuating magnetic field h^. Such a field would couple here to the magneti- 
zation Jl + J/j: integrating out would then induce a coupling — (Jl + JrY which gives 
us the Hamiltonian 

Hf^e = + 4) - AJl ■ Jr. (6.2) 
The model Hp_,E can also be used to describe the onset of ferromagnetism in an itinerant 



Luttinger liquid. It has in fact been examined earlier by Affleck where he obtained it 
as an effective model for the spin degrees of freedom in a = 1 Hubbard model. Affleck 
examined the RG flows of A for small A and obtained 

. In the mapping from the repulsive Hubbard model, and also from our rotor model, the 
initial sign of A is positive. From ( |6.3|) , Affleck concluded that A is irrelevant for all positive 
A, and that all such systems flow into the A = flxed point. We believe this conclusion is 
incorrect. It is clear from our arguments that for A large enough, Hp^E should undergo a 
phase transition to ferromagnetic ground state. This suggests that there is a critical value 
of A = Ac (with Ac > and of order unity) such that only systems with A < Ac flow into the 
A = flxed point. Systems with A > Ac are in the ferromagnetic phase E. The nature of 
the critical point at A = Ac, which controls the transition from phase F to E, is unknown: 
determining its structure remains an important open problem. 



VII. CONCLUSIONS 

This paper has introduced and analyzed the simplest model with Heisenberg symmetry 
which exhibits zero temperature phase transitions, and whose phases contain a net average 
magnetic moment. The model contained only bosonic quantum rotor degrees of freedom 
and offers the simplest realization of a quantum transition with an order parameter which 
is also a non-abelian, conserved charge. The analysis focussed primarily on the Hamiltonian 
H in ( |1.3|) , although variations were also considered. The results are summarized in the 
phase diagrams in Figs 

Some important properties of these phase diagrams deserve reiteration. Notice that for 
d > 1, there is no phase which is simply a non-quantized ferromagnet, with no other broken 
symmetry. Phase D has a non-quantized ferromagnetic moment, but it has an additional 
long-range order in the n fleld in a plane perpendicular to the ferromagnetic moment. We 
believe this is a generic feature of insulating spin systems: ferromagnetic ground states either 
have an integral or half-integral magnetic moment, or have an additional broken symmetry. 
Only m d = 1 does a phase like E appear: it has a non-quantized ferromagnetic moment, 
and is invariant under rotations about the ferromagnetic axis. However, even in d = 1 there 
is a remnant of the broken symmetry in the n fleld perpendicular to the moment: correlators 
of n have a power-law decay in space, and the linearly-dispersing gapless spin-wave mode 
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is still present (in addition to the usual quadratically dispersing ferromagnetic mode). It 
is also interesting to note that metallic, Fermi liquids of course have no trouble forming 
non-quantized ferromagnets; this is in keeping with popular wisdom that Fermi liquids are 
"effectively" one dimensional. 

A second interesting property of the phase diagram was pointed out in Section |I|: contin- 
uous zero temperature transitions in which there is an onset in the ferromagnetic moment 
only occur from phases which have gapless excitations. Thus there is such an onset from 
phase C (which has gapless spin waves) to phase D, but no continuous transition between 
phases A and D. 

We also examined the critical properties of the second order phase transitions in the 
model. In several cases, the critical theories turned out to be variations on the theme 
of a simpler quantum phase transition: the onset of density in a Bose gas with repulsive 
interactions as its chemical potential is moved through zero. This quantum transition had 
been studied earlier [^,^]. Because of its central importance, we obtained some additional 
results on its T = universal properties in Appendix 0. 

In the remainder of this section, we remark on issues related to those considered in this 
paper, but which we have not directly analyzed here. 



1. Itinerant ferromagnets 

We discuss implications of our results for quantum phase transitions in ferromagnetic 
Fermi liquids. As noted in Section |, there are two quantum transitions in this system, and 
we will discuss them separately: 

(z) Consider first the transition from a fully polarized Fermi liquid to a partially polarized 
Fermi liquid. This is rather like the transition from phase B to phase D in c? > 1, and 
the transition from phase B to phase E in = 1. Let us suppose that all electrons in the 
fully polarized state are polarized in the 'up' direction. Then an order parameter for the 
transition is simply the density of 'down' spin electrons. Along the lines of the analysis in 
the rotor model, we can derive an effective action for the down spin electrons simply by 
integrating out the up electrons. It turns out that the up electrons only mediate irrelevant 
interactions between the down electrons: the effective action for the down electrons is simply 
that of a dilute gas of free spinless fermions [0,0 . A possible four-point interaction like the 
\ip\'^ term for the Bose gas is prohibited by Fermi statistics; further a singular term, like that 
appearing in the models without V symmetry in Section [IV B l|b, appears to be prohibited 



here because the large Fermi momentum of the up electrons inhibits strong mixing with 
down electrons via emission of ferromagnetic spin waves. The properties of the free spinless 
fermion model are of course trivial, but it is quite useful to re-interpret them in the language 



of a quantum phase transition . It is also interesting to note here that the critical theory 



of free spinless fermions in = 1 is identical to that of dilute interacting bosons for the 



B to E transition m d = 1 p9[. This is in keeping with our assertion that in d = 1, the 
transitions in itinerant fermion systems are in the same universality class as those in certain 
rotor models. 

(m) Consider now the onset of ferromagnetism in an unpolarized Fermi liquid. A theory for 
this transition for d> 1 was proposed by Hertz 0] . This transition is similar to the transition 
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between phases C and D studied in Section |IV A| . In the end, our analysis used a method 
very close in spirit to that use by Hertz: simply integrate out all gapless modes not directly 
related to the order parameter. We have provided here some independent justification for 
such an approach in the rotor model, and our results provide support for the correctness of 
Hertz's analysis for d > 1. Precisely in d = 1, we have no theory for the transition between 
phases F and E, but we have argued that it should be in the same universality class as that of 
the onset of ferromagnetism in a Luttinger liquid of itinerant electrons. The critical theory 
of this transition is the main remaining open problem in the theory of phase transitions in 
quantum ferromagnets. 



2. Finite temperature 



Essentially all the analysis of this paper was at T = 0, and it would be interesting 
to extend it to finite T. The finite T properties of the quantum paramagnetic phase A, 
the Neel phase C, and the transition between them, have already been studied in some 
detail pT| , P^ . More recently, a field-theoretic analysis of the finite T properties above a 
quantized ferromagnetic state B has appeared |16[. It remains, therefore, to study the 
canted phase D and its phase transitions at finite T. We have shown in this paper that 
the transition between the B and D phase is similar (although not identical) to that in a 
dilute Bose gas [pSj^ — the finite T analysis should therefore be related to that in Ref |^ . 
The finite T properties of the phase D itself should also be quite interesting, especially in 
d = 2 where the long-range order disappears at any finite T. This is also the case of direct 
relevance to the quantum Hall effect experiments of Ref |TT]. The phase D could be modeled 
e.g. by the action ([4.25|) , supplemented by the additional constraint ri^ = 1, which can be 
imposed once we are well away from phase B. 



3. Effects of randomness 



At T = 0, but in the presence of randomness, one might expect that between the quan- 
tized ferromagnetic (B) and canted phases (D), there occurs a "canted glass" phase with a 
non-quantized magnetic moment, short-range mean Neel correlations, but a diverging mean 
Neel susceptibility; this can be seen by arguments analogous to those in Ref Further, 
general arguments imply that for arbitrarily weak randomness in c? < 2 and for suffi- 
ciently strong randomness in c? > 2, all the first order lines in our phase diagram are replaced 
by second-order transitions (implying, for instance, the absence of direct transitions between 
the quantized ferromagnetic phases). 
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APPENDIX A: MAPPING BETWEEN HEISENBERG SPIN AND ROTOR 

MODELS 



Consider a 'double layer' model of spin S Heisenberg spins Siaf^, Sib^ on the sites z of a 
d dimensional lattice; a, b are two 'layer' indices and = 1, 2, 3 are vector components. We 
study the Hamiltonian 

Hdi 

(Al) 

where < ij > is a sum over nearest-neighbor pairs. We will consider the case where G is 
antiferromagnetic (G > 0). Then, neglecting the inter-site terms, each site has a tower of 
states with total angular momentum £ = 0, 1, 2, . . . 2S. This tower is very similar to that of 
a single quantum rotor, the main difference being that the latter does not have an upper 
bound on its allowed angular momentum. For the low energy properties of interest in this 
paper, the upper bound is not expected to be important. Further, by constructing an on-site 
Hamiltonian which is a polynomial in (L^ is the rotor angular momentum), it is possible 
to mimic the actual eigenenergies of the tower of states in the Heisenberg system. 

Now consider the inter-site terms in H^i- For the on-site tower, the matrix elements of 
the rotor operators and are similar to those of Sa^ + Sb/j, and f{San — Sb/j,) respectively 
(/ is some constant); this correspondence becomes exact in the semiclassical theory. So 
we perform the replacement = (L^ -|- /n^)/2, Sb^ = (L^ — /n^)/2 in Hai- This yields 
precisely the inter-site terms in (|1.3|) with 

jy^ Ja Jb ~^ "^Jab _ j n2 Jo, ~^ Jb '^Jab _ -at ^Ja Jb f A2l 

Note that in systems with a layer interchange symmetry (J^ = J;,), the coupling M vanishes. 
Thus this interchange symmetry is equivalent to V. 



APPENDIX B: MEAN FIELD THEORY 
1. Rotors with gr = 

The Hamiltonian Hmf of Section || was diagonalized by determining its matrix elements 
in the basis of spherical harmonic states: |£, m) with £, m integers satisfying —i<m<i 
and i > 0. The matrix elements of the operators in this basis can be expressed in terms of 
Clebsch- Gordon co-efficients; the non-zero matrix elements are: 

(i, m\Lz\im) 
{i,m + l\L+\im) 

{£, m|n^|£ + 1, m) 



m 



((£ + m+ l)(£-m))i/2 
f {£ + m + 1){£ - m + 



(2£+ l)(2£ + 3) 



33 



{i,m+ l\h+\i+ l,m) 
{£, m — l\n_\i + 1, m) 

and their complex conjugates. Here, as usual, L+ = Lr^ + iLy and similarly for n^. All 
matrix elements not obtainable by complex conjugation of the above are zero, notice that 

only has non-zero matrix elements between states with angular momenta £ and £ ± 1. 
In particular, matrix elements of between states with the same value of £ vanish; this 
happens because is odd under P, and the |£, m) states have definite P-parity. 

The single-site Hilbert space has an infinite number of states, but in practice it was 
found that good accuracy was obtained by truncating the states above a maximum value of 
£ ^ 15. The results of the numerical calculation were discussed in Section | and in Fig|l|. 



{i-m){^t-m+l)\'- 
(2£+l)(2£ + 3) ) 

( {^ + m){^ + m + l)\ 
\ (2£+l)(2£ + 3) J 



2. Rotors with q = l/2 

The mean-field analysis proceeds in a manner similar to Sec We now have to introduce 
two distinct mean field Hamiltonians Himf, -f^2m/ for the two sublattices, each with their 
own effective fields A^i^, /ii^ and N2^, h2^. The expectation value of H in the ground state 
of the mean-field Hamiltonians is then minimized with respect to variations in these four 
effective fields. The numerical diagonalization of the mean- field Hamiltonians requires the 
matrix elements of the operators in the states of rotor with a monopole q. We determined 
these for g = 1/2 from Ref |32] and applications of the Wigner-Eckart theorem. The matrix 



elements of the are still given by those in ( pi| ) [i.e. they are independent of q) while 
those of the are 

{l/2,i,m'\h,,\l/2,im) = --^^-^^{1/2, i,m'\L,\l/2, im) 
(l/2,£,m|ni^|l/2,£+ l,m) 
(l/2,£,m + l|rii+|l/2,£+l,m) 



Ei f{i + m + l){i-m + l){2i+l)Y^'^ 



2£ + 2\ 2 y 

Si / {i - m){i - m + l){2i + 
2£ + 2 [ 2 



All matrix elements not obtainable from the above by complex conjugation are zero. Notice 
that, unlike (HI), now has non-zero matrix elements between states with the same value 



of i, and these matrix elements are proportional to those of L^; this happens because the 
\q > 0,i,m) states no longer have definite V parity. The results of this calculation were 
discussed in Section 0. 
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APPENDIX C: COMPUTATIONS FOR NEEL (C) TO CANTED (D) 

TRANSITION 



We will discuss the case without V symmetry, described by the action ([4.18|) in Sec- 
tion [iV A2| . The results with V symmetry in Section |1V A 1| follow as the special case 



Ao = 0. 

At the critical point, the propagator of ( [4 .181) is (dropping the subscripts on the cou- 
plings) 



r 



(6(A2 + 4))i/2(47r|r|)'^/2 \ xT ^-^ 
where 

-^i • (C2) 

The order v? contribution to the self energy is therefore 

u}) = 2u^ J d'^xdrG^x, r)G(-x, -r)e-^(^-™) (C3) 

After a lengthy, but straightforward, evaluations of the above integral we find, to leading 
order in e = 2 — c?: 

E(*=..)-E(0.0) = -g-5^i ({^^.jp^.^g/, + ... (C4) 

The renormalization constants Z, Zf^, and Zx in (|4.19| ) follow immediately from the above 
result. 

Similarly, the four-point vertex is 

u-u^ j d'^xdrG'^ix, r) - Au^ J d'^xdTG{x, t)G{-x, -r) 

+ 20 

~ Anb^/h (A2 + 4)3/2 + • • • ' lC5j 

which leads to the result for Z^ in ( [4.19|) . 

The computation of Z2 is very similar, and details are omitted. 

APPENDIX D: DILUTE BOSE GAS BELOW TWO DIMENSIONS 

In this Appendix we will study the Bose gas described by the action 

S = J d'xdr + + r|^r + . (Dl) 
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This action undergoes a T = quantum phase transition, at r = 0, which plays an important 
role in the models considered in the main part of this paper. In d < 2, this quantum 
transition obeys a no scale-factor universality |^9[ which we shall study here in greater 
detail. One consequence of this enhanced universality in c? < 2 is that the zero temperature 
density, n = of this Bose gas obeys 

n = C,0(-r)|r|'^/^ (D2) 



with Cd a universal number |2^. Here we shall show how to compute Cd in an expansion 
in powers of e = 2 — rf. It is known that the leading term is = 1/ [Aire) , and we shall 
explicitly determine the next term. We shall also compute the effective potential of the Bose 
gas to order e. 

We characterize the interactions in S by the bare dimensionless coupling qq defined by 

go = fi'^SdU (D3) 

where /i is a momentum scale, and recall that Sd = 27r'^/^/((27r)^r((i/2)). We define a 
renormalized dimensionless coupling g in a similar manner by replacing u with the value of 
the exact two-particle scattering amplitude at zero external frequencies and equal incoming 
momenta p, with = s^/i^. The dimensionless number s is arbitrary, and no universal 
quantity should depend upon it; we shall keep track of the s dependence as a check on 
the universality of our final results. The two-particle scattering amplitude is given by the 
sum of a series of ladder diagrams which can be evaluated exactly; this gives us an exact 
relationship between g and go 



^ = ^0 1 + 



gofi' f d^k 1 \ ' 



Sd J {27rYk^ + {k + 2py 



-1 



^ofl + ^l (D4) 



with 

^^^^im^^=l^m,^o(.-) (D6) 

We can now deduce the exact /5-function of g, 

Pg = -tg + Adg\ (D6) 

which has a fixed point aX g = g* = e/Ad- 

The effective potential (this is the generating functional of the one particle 

irreducible vertices) can be easily determined from the standard Bogoluibov theory of the 
ground state energy of a dilute Bose gas: 

rW =^\^\" + l\^\" + \j ^d [((^' + ^ + 2«IV^I')' - -{k' + r + 2«|^p)" 

(D7) 
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in the one-loop approximation. To express this result in a universal form, we have to evaluate 
the above integral, express u in terms of the renormalized coupling g using (P3| , D4 ) and, 
finally, set g at its fixed point value g = g*. First, we write ( p7| ) in the form 



2e 

Sd 



4e 



1 + ^? 



9 



.(k' + 2erf^'yigm+2) 



(P + 2er/x~7(^7H2) + 2)^ 



1/ 



+ 



Ad9' 



(D 



where ip = (2e/ SdY^'^'ip, and we expressed 5^0 only upto second order in g. The above integral 
can be evaluated in powers of e; as expected, the poles in e within the curly brackets cancel. 
Then, set g = g* = e/ Ad and expand the whole expression to order e. The /i and the s 
dependence disappears, and the resulting expression is completely universal. We can write 
it in the scaling form 



2e| 



-sgn(rj 



(D9) 



ySdlrl''/^^ J 

where $ is a universal scaling function, and all exponents are written in their expected exact 
form. Notice that there are no arbitrary scale factors in (|D9| ) — this is the no scale-factor 
universality of Ref [^. The scaling function $(?/) is determined by the above calculation 
to order e: 



7y' 



-sgn(y)(l + 2y){l + 4y + ?>y^Y'^ + l + Ay + ^\+ 



The condensate ^/'o = (V^) is determined by the condition dV ld'^\^=^f^ - 
total density of particles, n, is given by = dV /dr. Using ( P9| ) this gives us 



n 



sgn(r)|r|'^/2 ^ ^ ^^i^^^ yo^'{yo) = 



(DIO) 
0, while the 

(Dll) 



This gives a result for the density in the form (|D2|) with 



Cd — Sd — 
I 2e 



1 - log 2 



0{e) 



(D12) 



Finally, we recall that in c? = 1, the exact value of Cd is known 29 : Ci = Si. 



APPENDIX E: COMPUTATIONS FOR QUANTIZED FERROMAGNET (B) TO 

CANTED (D) TRANSITION 

At the critical point, the propagator of ( [4.38|) is (dropping the subscripts on the cou- 
plings) 
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G{x,t)= I 



{27tY 2tx -uj^M^X - iMouk^il + A) + k\l - b) 

j^(d-2)/2 I' g-Mox2/4A+r ^-Mqx'^ /4:X-\t\ 



-4+ Oir) - A_ -jj,^ ei-r) (El) 



where 



and 



. _ -((l + A)2-4A(l-6))V2±(i + A) 
2A 



We have assumed above, and in the remainder of this Appendix that A < 0. 
The order contribution to the self energy is therefore 

E{k, uj) = J d'^xdrG^x, t)G{-x, _r)e-^(^-^-'^") (E4) 

After a lengtliy, but straightforward, evaluations of the above integral we find, to leading 
order in e = 2 — d: 

E(fc, u) - S(0, 0) = -^^^(en.ib, A) - zMouR2{b, A)), (E5) 



and, as in ( P5|) , the renormalized four-point vertex 



u---——Rs{b,X), (E6) 
47rMo A e 



where 



Ri{b,X) 
R2{b,X) 



6A2((1 + A)2 + 3(1-&)(1-2A)) 
((1 + Xy - 4A(1 - 6))(2A2 + 9bX - 5A + 2)^ 

6A2(A + 36-2) 
(1 - 6)((1 + A)2 - 4A(1 - 6))(2A2 + 96A - 5A + 2) 
(1-Af -&(1 + 10A + 5A2) + 12&2A 
^> - (l-6)((l + A)2-4A(l-6))3/2 ■ ^^^^ 

To determine the renormalization of {ipl'^ insertions we need the vertex between a {ipl'^ oper- 
ator and a ip and a ?/'*; this is 

where 

= ((i + A)2-tAa-^'))^/^' ^^^^ 

These results lead immediately to the renormalization constants and (3 functions in Sec- 
tion lIVBlIb for A < 0. 
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FIGURES 




5 10 15 20 25 30 35 40 



ZJ 

FIG. 1. Mean- field phase diagram of H (Eqn |l.3D as a function of the couplings J and K at 
MZ = 4, ^ = 1, and a = 1/2; Z is the co-ordination number of the lattice. The model does not 
have V symmetry at these values, but the phase diagram for the V symmetric case is very similar. 
The mean-field result becomes exact in the limit of large spatial dimensionality, but the general 
features of the phase diagram are expected to hold for all d > 1. Thin lines represent second-order 
transitions while thick lines are first order. The quantized ferromagnetic phases have magnetic 
moment per site there is an infinite sequence of these phases for all integers £ > at larger values 
of and only the first two are shown. The phases have the following ground state expectation 
values, up to a global 0(3) rotation: (A) (L^) = 0, {n^) = 0; (B^) (L^) = (n^) / 0, (I/^.j/) = 0, 
{nx,y) = 0; (C) (L^) = 0, (n,) / 0, {n^,y) = 0; (D) (4,,) ^ 0, (n^,^) / 0, {Ly) = 0, {fiy) = 0. 
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FIG. 2. Mean-field phase diagram of with given by (|5.1| ), monopole charge q = 1/2, M = 
and other parameters and conventions are as in Fig|. Theph ases have the same physical properties 
as those of the corresponding phases in Fig|l]. The moment, of the quantized ferromagnet phases 

is now half-integral, and such phases exist for all half- integers I at larger K. The model is on a 
bi- partite lattice, and the expectation values (up to a global 0(3) rotation) in the various phases 
on the two sublattices, 1,2, are: (B^) (Li^) = {L2z) = i, (n^) = -{h2z) / 0, {Li^,iy,2x,2y) = 0, 

{■nix,ly,2x,2y) = 0; (C) (Li^) = -{L2x) / 0, (hix) = (mx) / 0, {Liy^iz^2y,2z) = 0, {niy^lz,2y,2z) = 0; 
(D) {Liz) = iUz) / 0, {fliz) = -{h2z) / 0, (Li^) = -{L2x) + 0, (ni^) = ljl2x) + 0, {Liy^2y) = 0, 

(n\y,2y) = 0; 
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FIG. 3. Expected phase diagram of W in d = 1 with monopole charge q = 0. The phases A and 
are as in Fig |l[ The d > 1 phase D becomes phase E in d = 1; the latter phase has (n^) = 0, 

{Lz) 7^ 0, {Lx,y) = 0. The magnetic moment is not quantized but varies continuously. The n^.^^ 

fields have algebraic correlations in space. 
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FIG. 4. Expected phase diagram oiTC in d = 1 with monopole charge an integer > 0. Phases 
are similar to those with the same labels in Fig ^ and ^. 
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FIG. 5. Expected phase diagram of W in d = 1 with monopole charge a half-integer q > 0. 
Phases are similar to those with the same labels in Fig ^ and ^. Phase F has no broken symmetry 
and algebraic correlations of all spin operators. Its low energy properties are described by the 1 + 1 
dimensional 0(3) non-linear sigma model with a topological term with co-efficient 6 = it. 
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